Unraveling Orbital Correlations via Magnetic Resonant Inelastic X-ray Scattering 
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Although orbital degrees of freedom are a factor of fundamental importance in strongly correlated 
transition metal compounds, orbital correlations and dynamics remain very difficult to access, in 
particular by neutron scattering. Via a direct calculation of scattering amplitudes we show that 
instead magnetic resonant inelastic x-ray scattering (RIXS) does reveal orbital correlations. In 
contrast to neutron scattering, the intensity of the magnetic excitations in RIXS depends very 
sensitively on the symmetry of the orbitals that spins occupy and on photon polarizations. We show 
in detail how this effect allows magnetic RIXS to distinguish between alternating orbital ordered 
and ferro-orbital (or orbital liquid) states. 

PACS numbers: 75.25.Dk, 75.30.Ds, 78.70.Ck, 74.70.Xa 



Introduction Ever since the seminal work of Kugel 
and Khomskii [l[ in the 1980s it has been known that 
orbital degrees of freedom can play a crucial role in corre- 
lated transition metal compounds. Orbital ordering and 
orbital-orbital interactions are not only closely tied to 
magnetic ordering and magnetic interactions, but orbital 
degrees of freedom have also been proposed to be of direct 
relevance to spectacular phenomena such as colossal mag- 
netoresistance in the manganites or superconductivity in 
the iron pnictides [2|-|4j. Yet, the precise nature of corre- 
lated orbital states, being either of ordered or liquid type, 
and their existence in different materials is intensely de- 
bated, which to a large part is due to the fact that orbital 
correlations turn out to be very difficult to detect exper- 
imentally. In fact, such experimental access would be of 
great help in unraveling the puzzling properties of many 
systems with orbital degrees of freedom, for instance the 
above mentioned iron-pnictide materials, where the type 
of the orbital ordering or its lack is heavily debated [B|-|7| 
or titanium and vanadium oxides where different theoret- 
ical scenarios — a rather exotic orbital liquid phase [8, 
or a classical alternating orbital-ordered state loL 11| - 
have been proposed. 

The experimental verification of orbital properties in 
correlated materials by neutron scattering is difficult be- 
cause neutrons are almost not sensitive to the orbital 
symmetries of the ground state, in particular in orbital 
systems the_ angular momentum is quenched by the crys- 
tal field 
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Traditional x-ray diffraction, instead, is 
dominated by scattering from the atomic core electrons 
while the resonant x-ray diffraction [l3, 14 1, particularly 



in the soft x-ray regime, the modern method of choice to 
detect orbital ordering, suffers from a very limited scat- 
tering phase space making Bragg scattering only possible 
for special orbital superstructures that have large spatial 
periodicities 15|. There being few orbital-ordering re- 
lated Bragg spots — if at all — leaves considerable room 
for controversies on the interpretation of experimental 
data 0-G1]. 

Recently, resonant inelastic x-ray scattering (RIXS) 



19h22| has been proven successful in measuring spin ex- 
citations in various cuprates I23rl28| . nickelates [2g|, and 
even iron-based compounds [30J . Here we show in a gen- 
eral setting how the polarization dependent intensity of 
magnetic RIXS directly provides an insight into the or- 
bital correlations in the ground state of correlated ma- 
terials. In particular, we verify that RIXS discriminates 
between different orbital states, e.g., the alternating or- 
bital (AO) order against the ferro-orbital (FO) order or 
the orbital liquid (OL) state. This method is applicable 
to any orbital-active material that has distinct disper- 
sive spectral features in its spin structure factor S(k, oj), 
for instance due to the presence of magnons arising from 
long-range magnetic ordering. 

RIXS cross section RIXS is particularly apt to probe 
the properties of strongly correlated electrons, for in- 
stance in transition metal (TM) oxides 22]. With an 
incoming x-ray of energy Wj n and momentum k; n an elec- 
tron is resonantly excited from a core level into the va- 
lence shell. At the TM L2.3 edges this involves a 2p — s- 3d 
dipole allowed transition. In this intermediate state, the 
spin of the 2p core hole is not conserved, as the very 
large spin-orbit interactions strongly couple the spin and 
orbital momentum of the core hole. A spin flip in the 
core allows the subsequent recombination of the core hole 
with a 3d electron that has a spin opposite to the electron 
that was originally excited into the 3d shell. The energy 
u; ou t and momentum k out of the outgoing x-ray result- 
ing from this recombination are then related to a spin 
excitation with energy uj = w ou t — Wj n and momentum 
k = k ou t k; n . 

The magnetic RIXS cross section at a TM £2,3 edge is 
in general 
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where e = e ln • (e out )'!' is the tensor that describes the in- 
coming and outgoing photon polarization, and H is the 
Hamiltonian describing 3d valence electrons with ground 
state |0) and energy Eq. The Fourier transformed tran- 
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FIG. 1. (color online) Schematic representation of the RIXS 
operator Oj^ e on a single site at the Cu 2+ Z/2,3 edge. To cal- 
culate the matrix elements of the operator between the same 
initial and final 3d orbital state, one needs to sum over all 
possible paths connecting them via a three step process, mul- 
tiplying at every step as indicated in the figure: (i) the in- 
coming polarization e a , (ii) — ici or the spin operator ±021^-7 
[positive (negative) sign for steps along (opposite to) the di- 
rection of the arrows], and (iii) the complex conjugate of the 
outgoing polarization e"jj (constants ci,2 depend on the edge). 



sition operator Ok, e = l/v^V^j Qj,e exp(ik • j) can be 
evaluated from the general expression for Oj e following 
the symmetry arguments in Ref. [HI [cf. Eqs. (8)-(10)] 



Oj,e 



'j<' 



Sj-We(dj), 



(2) 



D a j are the components of the dipole operator [22j , and 
Gj cx — 1C1+C2 Sj-IIj is the intermediate state propagator 
{c\ t 2 are constants depending on the resonant edge, see 
Fig.[IJ. The intermediate state transitions are expressed 
here by the operator ri 7j = J2 a p ^MP^ajPpj where £q/37 
is the Levi-Civita symbol and p* a ■ the creation operator 
of the 2p core hole in the p a orbital state. This com- 
pact expression for the core hole propagator leads to the 
schematic representation of the operator Oj e on a single 
site in Fig. [TJ 

While the intermediate state propagator Gj brings the 
spin dependence due to the spin-orbit coupling in the 
2p core hole states, the dipole operators D a - } act in a 
different way depending on the orbital occupancy on site 
j, so that the amplitude W e (dj) strongly depends on the 
orbital symmetry of the ground state at each site. Since 
this dependence is merely due to the properties of the 
dipole transitions and to the spin-orbit coupling, it is 
indeed generic to any TM £2,3 edge. 

Having analyzed the inherent dependence of the scat- 
tering amplitudes W e (dj) on the single site orbital oc- 
cupancy, we now investigate how the operator Oj e in 
Eq. ([2]) acts on the orbital ground state of the bulk. Here- 
after, we consider three different orbital ground states in 
a two-dimensional (2D) bipartite lattice (later we discuss 
a more general case): ferro-orbital (FO) order with the 
same a orbital occupied on each site, alternating orbital 
(AO) order with a (b) orbitals occupied on sublattice A 
(B), and orbital liquid (OL) ground state with the occu- 
pancies of a and b orbitals fluctuating similarly to the up 
and down spins in the spin liquid state. Thus we obtain 
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where Sj are spin operators, hjd are number operators 
for electrons in the 3d orbitals on site j, and where the 
vector amplitudes W e (dj) depend on the orbital symme- 
try dj of the ground state at site j. Here each component 
of the vector W e (dj) is a priori different and thus each 
spin operator is multiplied by a distinct amplitude, which 
can be related to the fundamental x-ray absorption cross 
sectionand therefore imp licitly depends on the orbital oc- 
cupancy dj at site j [2lL l3lL l32j . 

Orbital dependence of RIXS operator As stated above 
the orbital dependence of RIXS amplitudes W e (dj) is 
generic to any orbital system. Nevertheless, to be ex- 
plicit, we show how this dependence arises in the simple 
case of a Cu 2+ ion, i.e., with one hole in the Cu 3d or- 
bital. The amplitudes W e (dj) can be evaluated using 
Eq. ((2J as We(dj) oc (dj<T a |Oj ie |dj(T Q ) where \d- 3 a a ) is the 
state with a hole in the 3d orbital with spin a along the 
a axis. Since one needs here only to calculate the matrix 
elements of the operator Oj yf) on single site states, this 
can be done just by applying the dipole and fast collision 
approximations to the Kramers-Heisenberg formula for 

RIXS 33, 34j], so that Oj e = YlaB e ap^>\ jGjD a j, where a \ a i with aj being bosonic creation operators: keep- 



where the orbital pseudospin operator is T? = (hj a — 
1/2 for all sites j in the FO state 
i ±1/2 for every other site in the AO state, the 

operator 0- }>e acts differently on different orbital ground 
states. Below we show how this feature affects spectra, 
by calculating the cross section using Eqs. ([1]) and Q 
for six ground states with different orbital and magnetic 
configurations. 

FM systems with AO order We consider a 2D fer- 
romagnetic (FM) system with AO order (i.e., |0) 
I'M (gi AO)) with the spin interactions described by the 
effective Heisenberg Hamiltonian H = ^rSj with 

negative exchange constant J < 0. This spin-only Hamil- 
tonian follows from a Kugel-Khomskii spin-orbital model 
when the interactions between orbital degrees of freedom 
generating the AO ground state are integrated out (see 
Part 1 of the Supplemental Material). 

The spin wave (single magnon) excitation of such an 
ordered FM follows from the Holstein-Primakoff trans- 
formation for spins St = a 



aj and S? = 1/2 
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FIG. 2. (color online) Magnetic RIXS cross section J e (k, oj) 
for different magnetic (FM and AF) and orbital orders (FO 
and AO) along a high symmetry path in the Brillouin zone 
[where F = (0,0), X = (tt, 0), and M = averaged over 

incoming and outgoing polarizations. The FO (AO) order is 
formed by the x 2 — y 2 orbital (x 2 — z 2 and y 2 — z 2 ) while the 
spin quantization axis is in the xy plane. The color scale is 
nonlinear, since intensities of the AF spectra diverge at M, 
and at V in the AF-AO case. Spectra for the OL case (not 
shown) differ only quantitatively from the FO one. 



ing the quadratic terms in ay and Fourier transforming 
one obtains the bosonic Hamiltonian H = WkQ^ak 
with spin wave dispersion Wk = 2|J|(1 — 7k) where 
7k = (cosfc x -\- cos/by)/2. Furthermore one has Tj^jO) — 
exp(iQ • Rj)/2|0) where Q = (tt, tt) is the AO ordering 
vector, so that following Eq. Q one obtains 



AO)=±{[W e (a) + W e 



+ [W~ (a) - W- (b)] 4 +Q ||FM ® AO), (4) 

where W~ = Wg — iWg are the amplitudes for the spin 
flip transition, which can be calculated for the simple 
case of a Cu 2+ ion (cf. Fig. [IJ or for any other TM ion 



(cf. Refs.|2l|and|3J). Using Eq. (JD and the spin Hamil- 



tonian defined above, RIXS cross section can be directly 
calculated from Eq. (cf. Fig.f2]and Part 2 of the Sup- 
plemental Material). Due to the physical inequivalence 
of the two sublattices, the magnetic and orbital Brillouin 
zones are no longer the same, so that the backfolded 
branch of the magnon dispersion (pseudo optical branch 
in Fig. [5]) gains a finite intensity oc |W~(a) — W~(b)\ 2 
[cf. Part 2 of the Supplemental Material and Eq. (QJ] , as 
the spin flip amplitudes arc different for orbitals a and b. 

FM systems with FO order or OL state The above re- 
sult stays in contrast with the 2D FM case with FO order 
(|0) = |FM <g) FO)), for which one has f/|0) = 1/2)0) for 
all sites j. Again using Eq. ^ one obtains an equation 
for the operator Ok, e and for the cross section in Eq. fT}. 
In this case the orbital and magnetic Brillouin zones coin- 
cide since W~ (b) = W~ (a) and there is no pseudo optical 



magnon branch in the RIXS cross section (see Part 2 of 
the Supplemental Material and Fig. [2]). Finally for the 
2D FM case with an OL state (|0) = |FM ® OL)) with 
two fluctuating orbital states a and b the off-diagonal 
terms in Tj z |0) lead to orbital excitations and therefore 
can be omitted from Eq. ([3]), as we are interested only 
in pure spin excitations and not in coupled spin-orbital 
ones. Again the orbital and magnetic Brillouin zones are 
identical and only the acoustic branch is detectable (see 
Part 2 of the Supplemental Material). 

AF systems with AO order We consider a 2D antifer- 
romagnet (AF) with AO order (i.e., |0) = |AF ® AO)) 
with the effective Heisenberg interaction between spins 
as in the FM case but with J > 0. Similarly to the pre- 
vious case, the single magnon excitations are obtained 
by applying sequentially Holstein-Primakoff, Fourier, and 
Bogoliubov transformations and keeping only harmonic 
terms in bosonic operators a k and evk, cf. Ref. so 
that 



Ok.e | AF ® AO) = ~ { [W+ (a) + W~ (b)] u k <4 

~[W+(b) + W-(a)]v k al 

+ [W+(a)-W-(b)]u k+Q ai +Q 

+ [W+ (b) - W~ (a)] « k+ Qa k+ Q }|AF ® AO), (5) 

with W£ = Wg + iWg and where the Bogoliubov fac- 
tors are defined as u k — yv/2fik + 1/2 and vu = 
sgn(7k)\/ J/2^k — 1/2, and the AF spin wave dispersion 
is f?k = 2 J \J\ — 7 k . This form of the operator in general 
leads to a nonvanishing intensity when k — > V as a result 
of the AO ordering, see Part 2 of the Supplemental Ma- 
terial and Fig. In the case of ideal AF rik+Q = so 
that in contrast to the |FM ® AO) case one can observe 
only one branch in the RIXS spectrum (although any cor- 
rections to the Heisenberg model for which fik+Q ^ ^k 
will give rise to a pseudo optical branch in the spectrum, 
somewhat similar to the |FM ® AO) case). 

AF systems with FO order or OL state Again the 
above result stays in contrast with the 2D AF case with 
FO order (|0) = |AF®FO)) for which the RIXS operator 
has a simpler expression than Eq. ([5]) since W^(b) = 
W<f (a). In a similar way intensities for the 2D AF case 
with OL state (|0) = |AF ® OL)) are obtained (see Part 
2 of the Supplemental Material) . The intensity vanishes 
in both cases when k — >• T in agreement with Ref. 0, cf. 
Fig. [2] and Part 2 of the Supplemental Material. 

Discriminating different orbital states As shown 
above for FM and AF systems, RIXS spectra can dis- 
criminate an AO against FO order or OL ground states 
(cf. Fig. [2]). While in the FM case the pseudo opti- 
cal magnon branch signals the onset of the AO order, 
in the AF case the intensity of magnons with momenta 
k — > r does not vanish in the AO case, contrarily to 
the FO and OL case. This dependence is not due to dis- 
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tinct magnon dispersions for different orbital or electronic 
ground states [2, 3(| 37|, but to the orbital dependency 
of magnetic RIXS amplitudes. 

Furthermore, circular dichroism of magnetic RIXS in- 
tensities allows one to distinguish between different or- 
bital ground states, see Fig. [3] While for FM systems 
whether a circular dichroism is present depends on the 
symmetry of the orbital occupied, in the AF ones its pres- 
ence only depends on the system translational symmetry. 
Specifically, for |AF ® FO) (or |AF ® OL)) systems, cir- 
cular dichroism vanishes, while in the case of | AF ® AO) 
order (for which the RIXS spin flip amplitude is finite for 
both orbitals forming the AO ground state, cf. Ref. [20I ]) 
the circular dichroism is nonzero (Fig. [3]). 

In fact, if there is an AO order in a magnetic sys- 
tem, translational symmetry is broken into two phys- 
ically inequivalent sublattices. Consequently a pseudo 
optical branch in the magnon dispersion appears in the 
|FM ® AO) case. On the other hand, while a simple 
|AF Cg> FO) (or |AF ® OL)) system is symmetric under 
the combination of time reversal and a discrete trans- 
lation [38j], in the |AF ® AO) case the latter is broken. 
Macroscopically [39|, that means that the system is no 
longer symmetric under the combination of time rever- 
sal and translation. As a consequence, a finite circular 
dichroism appears, i.e., RIXS intensities (at fixed k and 
w) for left and right circular polarization of the incoming 
photon are no longer equivalent. 

Although the actual values of the W e (dj) transition 
amplitudes depend on the orbital symmetry at each site, 
differences in the RIXS spectra between the AO and the 
FO/OL ground states show up (cf. Fig. [2] and Fig. [3]) , as 
long as W e (a) 7^ W e (&). For this reason, the discrimi- 
nation between different orbital states does not rely on 
the particular orbital occupancy on the single site, but 
on the breaking of the translational symmetry caused by 
the onset of the AO orbital order. 

While other inelastic scattering methods have been 
theoretically proposed to detect orbital ordering 4oL 41 1. 
it should be stressed that, due to the onset of character- 
istic dispersion, the magnetic peaks in RIXS can, unlike, 
e.g., orbitons, be easily identified. Besides, as typically 
magnons interact weakly, quasiparticle peaks in RIXS 
spectra have sharp and well-defined line shapes which 
would rather not be obliterated by other low energy exci- 
tations (cf. Ref. [22j]) and their dependence on the orbital 
ground state is thus very pronounced. 

Conclusions We have shown in detail how ground 
state orbital correlations directly reflect themselves in 
magnetic resonant inelastic x-ray scattering (RIXS) in- 
tensities. It follows that measuring the RIXS spectra 
at transition metal £2,3 edges in correlated materials 
with orbital degrees of freedom and magnetic order, al- 
lows one to distinguish between different orbital ground 
states 42] . This is possible because in magnetic RIXS the 
spin flip mechanism involves a strong spin-orbit coupling 




FIG. 3. (color online) Circular dichroism D — {Ie L — 
I eR )/{Ie L + Ie R ) for RIXS spectra intensities at oj = f2k as a 
function of transferred momentum k for the AF state, where 
(br) is left (right) incoming circular polarization, for AO 
(FO and OL) state plotted with solid (dashed) line. 



deep in the electronic core so that, unlike in inelastic neu- 
tron scattering, the magnetic scattering spectra strongly 
depend on the symmetry of the orbitals where the spins 
are in. 

The method proposed here is of direct relevance to 
2D orbital systems, e.g., K 2 CuF 4 or Cs 2 AgF 4 with FM 
layers and predicted (but not yet explicitly verified) AO 
ordering 43l-l45|. as well as to three dimensional transi- 
tion metal oxides with orbital degrees of freedom such as 
LaMn0 3 , KCuF 3 , LaTi0 3 or LaV0 3 Q. In particular, 
in Part 3 of the Supplemental Material we predict mag- 
netic RIXS spectra for two different polytypes of KCuF 3 
with distinct orbitally ordered ground states. 
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SUPPLEMENTAL MATERIAL 

The following supplementary information consists of 
three parts: in Part 1 we discuss the origin of the spin- 
only Hamiltonians which describes the spin-orbital prob- 
lems in the main text of the paper, in Part 2 we give 
explicit equations for the magnetic RIXS cross sections 
as studied in the main text of the paper, and in Part 3 
we calculate magnetic RIXS cross section of KCUF3. 



1. Spin-only effective Hamiltonian 

In the main text we consider 2D systems with orbital 
order with the spin interactions described by the effective 
Heisenberg Hamiltonian 

(U> 

with exchange constant J. Such a spin-only Hamiltonian 
follows from a Kugel-Khomskii spin-orbital model [l| 
when the interactions between orbital degrees of freedom 
generating the AO ground state are integrated out (us- 
ing e.g., mean- field decoupling). In a typical AO and FO 
case with large crystal-field Jahn- Teller interactions this 
should always be possible 0. In the OL case this might 
be questionable but the experimental results suggest that 
also in that case it is a valid approach [3] . However, if this 
is not the case, then in principle it is not clear what is the 
nature of the elementary excitations in such systems [iJ . 
In fact, so far there is only one case,closely related to 
the recently studied problem in Ref. [4|, i.e., in the case 
of a FM-AO state with no Jahn- Teller interactions, in 
which the supposedly well-defined magnetic excitations 
dressed with the orbital ones should also give rise to two 
branches in the RIXS spectra. This is somewhat similar 
to the case described by Eq. in the main text but 
with the pure magnetic excitations replaced by the ones 
dressed by orbitons (which would mean that each of the 
two branches visible in RIXS in Fig. [2] in the main text 
would also be visible but have also some large incoherent 
spectrum and potentially different periodicity). 



2. Magnetic RIXS cross sections 

In Table U we give explicit equations for the magnetic 
RIXS cross sections for the three distinct ordered orbital 
ground states (AO, FO, and OL) and two distinct mag- 
netically ordered ground states (FM and AF). This is cal- 
culated by substituting Eq. (gj [or Eq. ©] into Eq. (JI} 
in the main text for the FM (AF) case respectively. Note 
that results presented in Table U are then used in the 
main text inter alia in Fig. [2] 



3. Magnetic RIXS cross section of KC11F3 

In what follows we calculate RIXS cross section for the 
magnetic excitations in the so-called A-AF state (i.e., FM 
planes coupled AF along the c direction) as stabilized in 
KCuF 3 below T < T N ~ 38K [1 S] . We consider here 
three different orbital states @: (i) C-AO state ( i.e., AO 
planes coupled FO along the c direction) as stabilized 
in the so-called (d)-type polytype of KCuF 3 below T < 
T s ~ 800K J3], (ii) G-AO (i.e., isotropic 3D AO state) 
as stabilized in the so-called (a)-type polytype of KCuF 3 
below T <T S ~ 800K 0, and (Hi) OL state (not realized 
in KCuF 3 but included for comparison). Furthermore 
we calculate the spectra for two different sets of {a, b} 
orbitals forming the above discussed ground states: (i) 
the 3d x 2_ y 2 and 3d 3z 2_ r 2 orbitals which can be preferred 
in the idealized case of vanishing interaction with the 
lattice (crystal field and Jahn- Teller interaction), and (ii) 
the 3d x 2_ z 2 and 3dj,2_ 2 2 which is the set of occupied 
orbitals preferred by the interaction with the lattice and 
which is probably close to the one realized in KCuF 3 0] ■ 

Similarly to the main text of the paper, we start with 
determining the spin wave excitations in the studied mag- 
netic structure. These originate from the Heisenberg-like 
Hamiltonian 



H = J X Y, Si-Sj + J 2 ^ S, Sj, (6) 

(UMM <y>iic 



where the spin exchange constants J\ < and J2 > 
lead to the onset of the A-AF ordered ground state. 
The anisotropic structure of this spin-only Hamiltonian 
stems from the full Kugel-Khomskii spin-orbital Hamil- 
tonian Q when orbital degrees of freedom (generating 
one of the above mentioned orbital ground states) are in- 
tegrated out. This means that the values of the spin ex- 
change constants J\ and J2 depend on the orbital ground 
state. 

As in main text of the paper, the spin wave excitations 
(single magnon) can be calculated by applying sequen- 
tially the Holstein-Primakoff, the Fourier, and the Bo- 
goliubov transformation and keeping only the harmonic 
terms in the bosonic operators and «k- Thus we ob- 
tain the bosonic Hamiltonian H = e^a^a-^ with spin 
wave dispersion = \J A\ — where A^ = 2|Ji|(l — 
7k) + Ji and i?k = «/2A*k with 7k = (cosfc x + cosfcj / )/2 
and /ik = cos k z . 

Moreover, using Eq. ([3]) in the main text we calcu- 
late the RIXS transition operator for a ground state cor- 
responding to a magnetic |S(Q)) state and an orbital 
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r M 


A T7 
Ar 


AO 


\W-{a) + W-{b)\ 2 8{uj-^) 
+ \W~(a) - W-(b)\ 2 5(uj -u; k+q ) 


|[W+(a) + W.-(6)]u k - [W+(6) + M/ e -(a)]« k | 2 5(a; - k ) 
+|[W+(a) - W e -(6)]« k+q + [W+(6) - W e -(a)H +Q | 2 5(w - fi k+ q) 


FO 


|Vy-(a)| 2 <S(^-o; k ) 


[|W+(a)| 2 + |W-(a)| 2 ]K - v^) 2 8(u - fl k ) 


OL 


|H/-(a) + M/ < r(6)| 2 <5(^-a; k ) 


[|W. + (a) + W+(6)| 2 + |W-(a) + V^-(fo)| 2 ](^ k - « k ) 2 S(aj - Q k ) 



TABLE I. Magnetic RIXS cross sections 7 e (k, w) for three different orbital (AO, FO and OL states) and two different magnetic 
(FM and AF) ground states, see main text and supplemental materials for further details. Constant factors are omitted. 



|0(Q)) one, with ordering vectors Q and Q respectively 

6k,e|S(Q)®0(Q)) = 

\{[W+(a)+W-(b) + W+(b)+W-(a)] 

x (u k - Uk)a k 
+[W+(a)+W-(a) - W+{b)-W-{b)\ 

x (u k+Q - «k+Q)a k+Q 
+[W-(a) + W-(b) - W+(a)-W+(b)] 

x K+Q + ^k+qK+Q 
+[W+(b) + W-(a) ~ W+(a)~W-(b)\ 

X K+Q+Q + ^k+Q+Q)"I+Q+Q}l S (Q) ® °(Q))' 

(7) 

where the Bogoliubov factors are defined as it k = 
v / A k /(2e k ) + 1/2 and w k = sgn(B k )^ k /(2e k ) - 1/2. 
Here the ordering vectors are defined through the equa- 
tions Sf |S(Q)> = exp(iQ • Rj)/2|S(Q)) and T/|0(Q)> = 
exp(zQ • Rj)/2|T(Q)) so that one has Q = (0,0, tt) for 
the A-AF phase and Q = (tt, tt, 0) [ Q = (tt, tt, tt) ] for 
the C-AO (G-AO) phase and Q = (0,0,0) for the OL 
phase. Note, however, that the above equations are valid 
also for other 3D spin and orbitally ordered phases with 
different ordering vectors and therefore can be used to 
study magnetic RIXS spectra for various other spin and 
orbitally ordered ground states. 

Finally, using Eq. ([7]) above and Eq. ([T]) in the main 
text, we can calculate the RIXS cross section, that reads 

7 e (k, lo) cc 

\W+(a) + W-(b) + W+(b) + W-(a)\ 2 

X (Uk - Vk) 2 S(0J - £ k ) 

+| W+ (a) + W~ (a) - W+ (b)-W-(b)\ 2 

x (u k+Q - v k+Q ) 2 5(uj - e k+Q ) 
+| W- (a) + W~ (b) - W+ (a)~W+(bf 

x (w k+ Q + w k+ q) 2 (5(a; - e k+Q ) 
+| W+ (b) + W- (a) - W+ (a) -W~(b)\ 2 

x ("k+Q+Q + u k+Q+q) 2 ^( w - e k+Q+q), (8) 



and is shown in Fig. for the above mentioned three 
different orbital ground states and for two choices of the 
orbital sets forming these ground states. 



While no clear signatures in the RIXS spectra allow 
one to distinguish between the two sets of single site or- 
bital occupancies {x 2 — y 2 , 3z 2 — r 2 } and {x 2 ~ z 2 , y 2 — z 2 } 
(a subtle intensity shift between the optical and the 
acoustic branch is hardly visible), differences between 
the C-AO, G-AO orders and OL state show off strikingly 
(see Fig.@|. In fact, in the OL spectra only the acoustic 
branch is present, whereas in the AO orders (G-AO and 
G-AO) spectra the pseudo-optical branch is present as 
well (compare differences between the FO and AO states 
in the main text). Moreover, the optical branch spectral 
intensities at M = (tt, tt, 0) discriminate between the two 
different AO orders: while they vanish in the G-AO, they 
diverge in the the G-AO case. 



Therefore, magnetic RIXS cross section allows one 
to distinguish between various orbitally ordered phases 
which are predicted to be stable in KCUF3 [3]. Further- 
more, the spectrum of the OL phase (which is not stable 
in the magnetically ordered phase of KCUF3 Q) is strik- 
ingly different than the one of the ordered phases which 
shows that the magnetic RIXS cross section strongly de- 
pends on the orbital correlation - as already discussed in 
detail in the main part of the paper. 
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FIG. 4. (color online) Magnetic RIXS cross section J e (k,oj) for KCUF3 along a high symmetry path in Brillouin zone for three 
different orbital ground states (from left to right, C-AO, G-AO and OL state), formed by alternating x 2 — y 2 and 3z 2 — r 2 
orbitals (top row) , and alternating x 2 — z 2 and y 2 — z 2 orbitals (bottom row) , and with the spin quantization axis in the xy plane. 
The cross section is averaged over incoming and outgoing polarizations and we assume that |«7i | = 0.06 in KCUF3 [3, [(|. 
The high symmetry points in the Brillouin zone are defined as Y = (0, 0, 0), X = (n, 0, 0), M = (ir, n, 0), and Z = (0, 0, it). The 
color scale is nonlinear, since intensities diverge at Z in every case, and at M in the G-AO case. 



